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We introduce a sequential model for the deposition and aggregation of particles in the 
submonolayer regime. Once a particle has been randomly deposited on the substrate, it 
sticks to the closest atom or island within a distance I, otherwise it sticks to the deposition 
site. We study this model both numerically and analytically in one dimension. A clear 
comprehension of its statistical properties is provided, thanks to capture equations and to 
the analysis of the island-island distance distribution. 

§1. Introduction 

The morphology of a system growing by deposition and aggregation of particles^) 
depends on several factors and it is not exaggerated to say that a detailed description 
would require a theory, or at least a model, for each different physical system. In 
spite of this, the use of simple and general models for studying growth processes is 
widespread. This is justified by two big reasons: First, such models allow to study 
large-scale features which are common to different physical systems; Second, they 
allow to connect apparently different phenomena. 

In the following we are referring to the growth process of a crystal surface by 
molecular beam epitaxy,^) with a special interest in the submonolayer regime,^) 
where only a fraction of a monolayer has been deposited. With these caveats in 
mind, a simple and very popular model for the growth process^-* accounts for random 
deposition of atoms and their thermally activated diffusion till they meet irreversibly 
another atom (nucleation) or island (aggregation). The system is therefore made up 
of diffusing particles (adatoms) and still islands. New atoms, which are continuously 
provided by deposition, can not leave the surface, because evaporation is forbidden 
at low temperature. Hereafter, this model will be called 'full diffusion' model. 

The diffusion length, ij^, is an important quantity. It measures the typical 
(linear) distance walked by an adatom before being incorporated.*^^ It depends on 
the deposition rate, F, and the hopping rate, h, through the relation £d ~ (/i/F)^/'*', 
with 7 = 2{d* + 1), where d* is the spatial dimension of the system. Once an atom 
has been deposited it visits a region of linear size and sticks inside that region 
after a time ~ ^u/h. The probability that another deposition occurs during this 
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time in that region of size is 

FtX' - {F/h)il+^' ~ el^''*-^ ~ i-^* . (M) 

Since £d ^ 1 such probabihty is very small, meaning that the sticking process 
of a new adatom is not influenced by atoms deposited later on. These arguments 
lead us to consider an even simplified model, which we are now going to describe. 

Once an atom has been deposited, it immediately searches the surroundings for 
another atom or an island. If it exists, it sticks to it, otherwise it sticks to the 
deposition site. Next, another atom is deposited. The search will be implemented 
deterministically through a capture area: the newly deposited atom (adatom in the 
following) looks for the closest atom/island within a distance i. Each adatom will 
therefore attach to an existing atom/island or it will stick to the deposition site and 
collect atoms deposited later on in its capture area. It is to be noted that in the 
present model an adatom moves around only at the time of deposition, but never 
afterwards. It is also important to stress that the size of the effective capture area 
reduces in the course of time, if newly deposited atoms have their own capture area 
overlapping with it. The reason is that particles at a distance d, with ^ < d < 2^, do 
not capture each other, but their capture areas overlap, and their capture distances 
along the joining line are reduced to d/2. 

We are going to study the above model, both numerically and analytically, in 
one spatial dimension {d* = 1). This choice is due to the possibility to provide a 
full theoretical comprehension of numerical results, using two main analytic tools: 
capture equations and the analysis of the island-island distance distribution. 

In the literature of submonolayer deposition, two models for representing islands 
are usually found: 'extended-islands' and 'point-islands'. In the former case, the 
physical size of an island increases proportionally to the number of adatoms attached 
to it and its shape depends on additional factors. In the latter case, islands have a 
physical size equal to a lattice site whatever is the number of atoms they are made 
up. The two models are statistically equivalent at small coverage (see next Section). 
We will use a 'point-island' model because of its simplicity. 

A last remark relates to the physical interpretation of the capture length and to 
the novelty of our model. If we refer to a system which displays a thermally activated 
diffusion process, our capture length can be read as an effective way to implement it. 
The idea is not new: it has been used, e.g., by Michael Biehl, Wolfgang Kinzel and 
coworkers in several papers^) for studying multilayer growth and it is related to the 
model of random deposition with surface relaxation.^) However, we are not aware 
of the application of the idea of capture length to study the submonolayer regime. 
On the other hand, the capture length may represent nonthermal post deposition 
processes, as found, e.g., during the adsorption of rare gas atoms on metal surfaces. ''^ 

This paper has three Sections, beyond the present Introduction. In the next 
Section we provide a summary of the main numerical results, while in Section |21 we 
derive the analytical results and compare them to numerics. In the final Section we 
discuss the hypotheses underlying our model and suggest a possible interpolation 
between our model and the 'full diffusion' model. 

Some remarks on notations are in order. In the analytical calculations we gen- 
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erally use 'dimensional' quantities, for example for the distance d between islands. 
However, when we display numerical results and compare them with analytics it may 
be useful to use reduce variables, for example x = d/i. In order to avoid any pos- 
sible misunderstanding, we say that d,d,doo,dnn,dnnn,di,dr,y are all dimensional 
distances. In particular, d is the average distance between neighbouring islands and 
doo is its asymptotic value for large coverage (i.e., for large deposition time). The 
only adimensional distances that are used throughout the paper are x and d = doo/i- 
Finally, two probability distributions for the distance will be introduced. Pi 
and P2. They refer to the distance between first (nn) and second (nnn) nearest 
neighbouring islands, respectively. 

§2. Numerical Simulation 

Our model is defined as follows. We choose randomly a lattice site and check if 
some atoms or islands exist within a distance i. If they exist, the newly deposited 
atom sticks to the closest atom/island, otherwise it sticks permanently to the depo- 
sition site. When an atom sticks to another atom, we get a new island. 

The model has two main features: it is a sequential model and random deposition 
is the only source of noise. If Na atoms are deposited on a substrate of L sites, the 
ratio 9 = Na/L defines the coverage 9, which is also equal to the product of the 
deposition rate and the deposition time, 9 = Ft. The statistical properties of the 
model depend on the product between the capture area c = 2l+\ and 9: this product 
is the average number of particles deposited in the area c. All simulations have been 
done for £ > 10, so that in the following we can write c ~ 2£ and introduce the 
parameter p = 9i. Point-island and extended-island models have the same properties 
at small coverage, let's say 9 < 0.2. This condition can be satisfied at any p, if the 
capture length is large enough, £ > 5p. 

Let us now give a qualitative description of the growth morphology followed by 
a summary of numerical results. Since all quantities depend on p = 9i only, we can 
think to keep £ fixed and vary 9 or the other way round. For pedagogical reasons we 
are describing the growth morphology with increasing 9. 

At the very beginning, deposited atoms do not interact and each atom or island 
grows in size according to the number of particles deposited in its capture area 
c = 21. The average distance between atoms or islands is very large. With increasing 
9, atoms or islands densities increase and their capture areas start to overlap. From 
then on the growth processes of different islands are no more independent. The 
average distance is now between i and 2£. At large 9, almost all neighbouring islands 
have distances d < 2£. The majority of islands therefore increases in size according 
to the number of particles deposited in a capture area equal to half the distance with 
the left and right neighbours. Sometimes two neighbouring islands have a distance 
d > 21 and the interval may be empty or contain one atom. Densities of atoms and 
islands change because empty intervals are filled by atoms and atoms are transformed 
into islands. Let us now sum up the numerical results with the aid of the figures. 

Fig.^summarizes all the main average quantities: the densities of atoms (circles) 
and islands (stars) in the main figure and the average size of islands in the inset. At 
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Fig. 1. Numerical (symbols) and analytical (lines) results for the adatom (circles) and island (stars) 
densities. At small p (full lines) adatom density varies as fii = pexp( — 2p) and island density 
varies as fiis = 5 — (p + 5) exp(— 2p). At large p (dashed lines), fii decays as and and fiis 
approaches the asymptotic value 1/p (see Eas. l3-27l3-28ll '). Inset: The average size of 

islands as a function of p. Full line: small p approximation. Dashed line: large p approximation 
{sav = dp). All analytical results do not contain any fitting parameter. 

small p, the density of atoms (hi) increases linearly, of course, and the density of 
islands {fiis) increases quadratically (the tilde over a density means a reduced density, 
np = £n[], see next Section). This behaviour is related to the specific features of the 
nucleation process in the present model, which is due to the deposition of an atom in 
the capture area of an existing atom and it is not mediated by surface diffusion, as in 
epitaxial growth. In the latter case his increases with 6 (i.e., withp) as his ~ ^"^i with 
an exponent*) x = 3 in two dimensions. Mean- field rate equations predict the same 
exponent in one dimension as well, but in this case the assumption dhis/dO « h\ is 
wrong. The correct theory for nucleation on top of a terrace^-* suggests X = §• 

At large p, his saturates but does not decrease because coalescence of islands is 
not possible in a point-island model. The limiting value, n^, is related to the reduced 
asymptotic average distance between islands, d = d^ol^-, by the trivial relation = 
Xjd. Island density converges to h°^ as while adatom density decays as 
These behaviours (see next Section) are due to the existence of intervals slightly 
larger than 2£, which are 'filled' by an adatom, which is subsequently transformed 
in a island. The probabilities of the two processes differ: filling an interval d > 2£ 
with an adatom requires deposition on a small region of size d — 2i, while deposition 
on a region of size d^o is enough to make an island from an adatom. The process 
(adatom) — > (island) is therefore quicker than the creation of new adatoms: this is 
the reason why hi decays more rapidly than his converges to h^. A quantitative, 
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more rigorous analysis can be found in Section f3.3l 

The inset of Fig. Ogives the average size of islands, Sav It is of order two for 
small p. At large p, deposited atoms are shared ciniong existing islands, and. Sav 
increases linearly with 0, i.e. with p. 
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Fig. 2. The size distribution of islands for different values of p (Main: p = 0.2, 0.5, 1; Inset: p = 2, 5). 
Symbols: numerical results. Lines: analytical results according to the small p approximation, 
Eq. ll^ . 



Fig. 12 and Fig. |31 plot the size distribution of islands, at small (Fig.Ej) and large 
(Fig. IS)) p, respectively. At small p, when different islands are almost independent, 
the size distribution is expected to follow a Poisson distribution: this is confirmed by 
lines, which reproduce well numerical data (symbols) till p < 2. At large p an island 
grows in size because it collects atoms deposited in a region of size ^{di + dr), where 
di^r are the distances to the nearest left and right neighbour. In this regime, the size 
distribution is therefore strictly related to the distance distribution between second 
nearest neighbouring (nnn) islands: this is proved in Fig. [21 by comparing the two 
(rescaled) distributions, plotted as full line and dashed line, respectively. The dotted 
line is the analytical nnn distance distribution {P2), as derived from the island-island 
distance ditribution (Pi) under the assumption that there is no correlation between 
neighbouring distances (see the discussion in the next Section). 

The previous discussion shows that the distance distribution plays an important 
role to understand the statistical properties of the model. In Fig. 0^ we report the 
distribution Pi{x) of the normalized nn distance x = d/i, for several values of p (note 
the log-scale on the y axis). Distances smaller than one (i.e., d < i) are forbidden, of 
course. Two regimes are clearly visible, for x smaller and larger than 2. The reason 
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Fig. 3. Numerical results for the size distribution (full line), as a function of 2s /p. Distance dis- 
tribution between next-nearest-neighbouring islands: numerical (dashed line) and analytical 
(dotted line) results {x = d„nn/i)- 



is straightforward: the capture areas of two neighbouring islands at distance x < 2 
overlap, while they do not if x > 2. This implies that new islands can nucleate in 
between only in the latter case. 

For X > 2, Pi{x) decays exponentially, Pi{x) ~ exp(— ax) with a prefactor a 
depending on p (see Fig. El- At smaller distances, x < 2, Pi{x) depends algebrically 
on X. Furthermore, it is an increasing function at small p and a decreasing one at 
large p. For very large p, it converges to a limiting shape (see Fig. E)), except for x 
very close to one. 

§3. Quantitative Analysis 

3.1. Densities and size distributions for small p 

In the following ni is the density of adatoms (adatoms per lattice site), Ug (s > 2) 
is the density of islands of size s, and His = '^s>2 total density of islands. 

The quantities tt-q = ^n^ are the 'reduced' densities and they mean the number of 
adatoms/islands per capture length. In the limit of small p deposited particles do 
not interact because their distance is larger than 21 and each adatom or island has 
a capture area equal to 2i. In this limit it is possible to write down the following 
capture equations (p = 9i): 

^ = 1 - 2{nis + 2ni) (3-1) 
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Fig. 4. Simulation (a) and theory (b) for the island-island distance distribution, as a function of 
X = d/l, for p = 0.2, 1, 10, 100. 



^ = 2ni (3-2) 

dp 

which can be easily solved, giving 

ni = pexp(-2p) , his = I - {p + I) eyip{-2p) . (3-3) 

These expressions, that are consistent for any p, are reported in Fig. ^ as full 
lines and compared with numerical results (symbols). Comparison is almost perfect 
till J? ~ 1. For larger values overlapping of capture areas is relevant and the above 
approximation is no more valid. The main defect of (|3-3|) . due to the assumption of 
a constant capture area, is that hi and 2 — f^is — fif^ — 

decrease exponentially 

for p 00, instead of that in a power law. The limit of large p will be treated in 
Section 13.31 

The equations for the density evolutions of islands of any size s > 2 are 

^ = 2{hs.i-hs) (3-4) 
dp 

and they can be solved recursively, 

n, = -^exp(-2p) s>l (3-5) 
2 SI 

showing that at small coverage the size distribution is Poisson-like. In Fig. |21 we 
compare the above analytical expressions with numerical results: they match very 
well till p « 1. The average size Sav of islands {s > 2) is calculated as 

„ _ET=2sns_^^ exp(2p)-l ^^^^^ 



his exp{2p) - (1 + 2p) ' 
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This small p approximation for Sav is reported in the inset of Fig. ^ as a full 
line. It also reproduces the correct linear behaviour at large p {sav = ap), but with 
a wrong factor a. The reason is that the small p approximation gives an asymptotic 
value for the island density, = ^ , which is smaller than the actual value (compare 
stars and full line in Fig. therefore a = is overestimated. 

Finally, expression (|3-5|) for hg might be slightly improved for p > 1 by replacing 
{2p) with a parameter A in Eqs. ()3-5() and H3-6p . This parameter should be determined 
by Eq. Q 

using the actual numerical value for Sa^. 

3.2. Island-Island distance distribution 

In Fig. 0^ we plot the distribution of distances between nearest neighboring 
islands as a function of the 'reduced' distance x = d/i. The existence of two regimes 
separated at x = 2 is due to the possibility of creating a new island within an interval, 
only if d > 2£. The weight of the two regimes changes with p, because the fraction 
of distances d < 2i increases from zero to one as p increases. 

For X > 2 the distribution has an exponential behaviour, Pi{x) ~ exp(— qx). In 
this regime the capture areas of two neighbouring islands do not overlap and Pi{x) 
is just the probability that no island nucleation occurs in time 6 in the interval of 
size d between the two islands. ^''^ If Lo{t) is the nucleation probability per unit time 
and unit length, at time t, neglecting correlations of nucleation events we obtain 



Pi(x) = exp 



d dtuj{t) 







(3-7) 
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Fig. 6. Numerical (circles) and analytical (full line) results for the 'asymptotic' (p — 20) distance 
distribution Pi{x) in the interval 1 < a; < 2, i.e. £ < d < 2i. 



Nucleation occurs when an atom is deposited in the capture area of an already 
existing atom. If a deposition event increases the (total) capture area by A, then 
duj/dt = (1 — Lo)A, i.e. = 1 — exp(— Z\t). Inserting the last expression in 1)3 -71) 
and introducing 6 = A/£, we obtain 

a{p) = r dT{l-e-^^) (3-8) 
Jo 

= p-{l-e-^P)/6. (3-9) 

In fact, the quantity 5 is not constant, but it decreases in time, i.e. with p, being 
equal to two for p — > 0. In Fig. [21 we compare numerical data for a with Eq. (|3-9|) . 
assuming 6 = 1 + exp(— p), which gives the correct limits 6^2 for small p and 5 ~ 1 
for large p.^^^ The agreement is fairly good. 

If X < 2 the determination of Pi{x) is more difficult because of the 'interaction' 
between the capture zones of the two neighbouring islands. The existence of two 
islands in y = and y = d after a time 6 can be depicted as follows. If we assume 
that the (first) atom in zero is deposited at time ti and the atom in d is deposited 
at time t2 > ti, we require that no deposition takes place in a region {d + 2i) during 
time ti, no deposition in a region d between ti and t2, at least one deposition in 2i 
between ti and 6, at least one deposition in d between t2 and 0}'^^ Integrating over 
ti and t2 the product of all the above probabilities provides Pi{d) in the hypothesis 
that neighbouring intervals have no influence, 

t-e re 
Pi{d) = / (itiexp(-2£ti) / dt2ex.p{-dt2) X 
Jo Jti 
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X [1 - exp{-d{e - ta))] [1 - exp(-2£(e - h))] 
= exp{-{d + 2i)e)[h/d- h] (3-10) 

where 

exp((2£ + d)e) - 1 exp(2£g) - 1 exp(M)-l 

= WT~d 2i d + ^ ^^-^^^ 

h = {2£9exp{2ie) - exp{2e0) - 2(^6^ + l) /if . (3-12) 

Prefactors are missing in Eqs. (|3-7l3-lUp . They are determined imposing the 
continuity of Pi{d) in d = 2i and its normalization in oo). The result of this 
procedure is shown in Fig. whose qualitative comparison with simulation results 
(Fig. is rather satisfying. The quantitative comparison is good for small p, but 
is not for large p. The reason is simple: for large p, the average distance d between 
islands is smaller than 2i, which means that neglecting the left and right neighbours 
of the two islands located in y = and y = d is no more correct. It is possible to 
take them into account in a 'mean field approximation', by assuming the existence 
of such neighbours at a fixed distance d. At small p, d > 21 and this refinement is 
inessential, but for large p is relevant. 

The calculation proceeds along the same lines leading to ()3-10|) with the differ- 
ence that now capture areas are limited by the presence of two additional islands in 
y = —d and y = d + d. If we limit ourselves to the case p ^ 1, the refined procedure 
gives 

^i(^) = + doc - 2e)id + 2doo - 2i) ' ^^'^^^ 
where Pq and doo should be determined self-consistently through the conditions 

r2£ r2e 



f 

Jl 



dyPi{y) = l , / dyyPi{y) = d^. (3-14) 



In Fig.iniwe compare the analytical and numerical results for Pi{d) in the regime 
d < 2i, using the reduced distance x = d/L A couple of remarks are in order. First, 
the analytical curve has no fitting parameter. Second, the inverse of the reduced 
average distance d = doo/^ has the value 1/d = 0.738, determined self-consistently 
from the above procedure. This value agrees very well with the asymptotic island 
density for large p, ~ 0.75. The value of d also allows to draw the dashed line 
in the inset of Fig. ^ which gives the theoretical prediction Sav = dp for the average 
size of islands at large p. 

A final comment concerns the shape of Pi{x) for x < 2, in the limit p — > oo. The 
curve with circles in Fig. El refers to p = 20. Numerical results show that a limiting 
shape does exist for x > 1.1, while Pi{l) seems to have a logarithmic divergence. 
This behaviour is due to distance-distance correlations which are not taken into 
account by the mean- field approximation leading to Eq. H3-13() . 

3.3. Size distributions and adatom/island densities for large p 

Fig. 131 shows that the asymptotic size distribution hg (full line) agrees well with 
the distance distribution P2{x) (dashed line) between next-nearest-neighbouring is- 
lands (x = dnnn/^)- The reason of that agreement is easily explained, because the 
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density of islands is almost constant for p ^ 1. Nearly all deposited atoms are 
captured by preexisting islands, which grow according to their capture area. The 
capture area of each island is just half the distance with its left neighbour plus half 
the distance with its right neighbour, ^{di + dr) = \dnnn- Therefore, for large p the 
size s grows accordingly to the relation s = dnnrfi 1"^ and the distance distribution 
between nnn islands is equivalent to the size distribution of the islands. 

In order to compare hs and ^2(2;), the size distribution is plotted as a function 
of Isjp. The two curves slightly differ for x < 2(s < p\. P2{x) vanishes for x < 2, 
while hs has a tail at small size. This tail disappears in the limit p — > 00, but very 
weakly, because his converges to as 1/p only (see Fig. P). 

It is interesting to compare ^2(2;) as derived from simulations with the nnn 
distance distribution, as derived from the nn distance distribution, Pi{x), in the 
hypothesis that neighbouring intervals are independent (Fig. |51 dotted line). From 
the general relation 

P2{d) = J dyPi{y)Pi{d - y) , (3-15) 



for large p we have 

P2{d) = pi j 



dy 



{y + doo- 2£){y + 2d^ - 2l){d - y + d^ - 2t){d - y + 2d^ - 21) 

(3-16) 

where Pq is the (known) normalization factor for Pi{y)- The qualitative behaviour 
of P2{d) can be understood as follows. For large p, Pi{y) is practically non- vanishing 
only in the region i < y < 2i, so that P2{d) does not vanish for 21 < d < 41. 
The shape of P2{d) is determined by two factors: Pi{y) is a continuously decreasing 
function, and different d have a different 'weight' p{d). If we rewrite H3-15|) as a 
two-dimensional integral, P2{d) = dyidy2Pi{yi)Pi{y2)S{d — yi — y2), the weight 

is just the quantity p{d) = J^^ dyidy2S{d—yi — y2), which has a symmetric maximum 
in d = 3£ and vanishes at the extremities d = 21, 4i. The non-analytic maximum is 
the responsible for the change of slope of P2{d) for d/l = 3, while p{2t) = p{M) = 
explain the vanishing of P2{d) at the same points. 

As d increases from two to four the functions Pi in (|3T5j) are evaluated, in 
average, at increasing values of the argument. Therefore, for d > 3i both the weight 
p(d) and the product of Pi are decreasing functions of d: this explains the fast 
decreasing of P2{d) in that interval. For d < 3i the weight is an increasing function 
vanishing in 2£ and the product of Pi is a decreasing function: this justifies the 
presence of a maximum of P2{d) in that interval. 

Comparison of dashed and dotted lines in Fig. |31 shows a major difference in the 
region dnnn ^ 2i, i.e. for dnn <L £■ This disagreement is due to two reasons. First, 
the theoretical Pi{d) underestimates the true island- island distance distribution for 
d c:^ i (see Fig. inj. Second, Eq. (|3-15p assumes there are no correlations between 
neighbouring intervals, which is not the case. 

Let us finally discuss the adatom and island densities in the limit of large p 
(Fig. n dashed lines). In the large p regime, the surface is a sequence of islands 
separated by distances d < 2i (with average value doo)- Rarely are there intervals 
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with d > 11: some of them, equal in number to /q; are void, the others (/i) contain 
one atom. The total number of intervals is approximately equal to It = L/doo and 
it is assumed to be constant, because It ^ Iq) ^i- The number of atoms is equal to 
h. 

If Na is the number of deposited atoms {Na = OL), Iq and Ii satisfy the following 
equations 

dio Ay 



dNa 

dh 



L 



-In 



Ay 



dr.. 



(3-17) 
(3-18) 



dNa L L 

where Ay is the 'active' region of an interval d larger than 2i (active means that 
a deposition event in such region creates a new adatom). We can evaluate it as 
Ay = {d)y — 21, where the average (d)> is performed on intervals d > 2£ only. 
For d> 21 (and large p) the distance distribution between islands is 

P^{d) = Pi(2£)exp[-^(d- 2^)] , (3-19) 

with Pi{2i) which can be determined by Eq. H3-13() . 

The integration of the previous equation gives the probability that a distance d 
is larger than 2£, 

P{d > 2e) 



Co 
P 



withco = [2dln(^)]-i 



2d-2' 



Much in the same way we can determine Ay, 



/2T dye- 

which gives Ay = 1/9. 

Eqs. (|3T7I3-18|) can now be rewritten as 

d/o Iq 



(3-20) 



(3-21) 



dp 
dh 
dp 



P 



P 



dh 



whose solutions are 



h{p) 
hip) 



Cl 

p 

cid 



(dp)' 



n=3 



in 



1)1 



{dpy 



cil_ 
d P^ 



(3-22) 
(3-23) 

(3-24) 
(3-25) 



Eq. H3-25() proves that the adatom density vanishes as 1/p^. Since Iq ~ 
island density and the total density (adatoms+islands) both converge to the asymp- 
totic value with corrections of order Finally, we can determine analytically ci, 
because 



Co 
P 



P{d > 2i) 



lo + h 
It 



cidpo 1 
L p 



(3-26) 
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so that ci = {L/doo)co and 

hi=£^ = ^ . (3-27) 

L 2d3ln(|g)P^ 

In Fig. ^ we compare numerical data (circles) with the previous expression for hi 
(decreasing dashed line): the slope (-2) is pretty correct, but the analytical prefactor 
is of order 0.23, to be compared with the numerical value, 0.34. Along the same 
lines, it is possible to determine his for large p: 

nis = -~ , , , - 1 + ^ . 3-28 

d 2dHn '^)P\ dp J 

\ 2d— 2 J 

This analytical expression is compared succesfully to numerical data in Fig. ^ 
(see the dashed line superposing to stars). 



§4. Comments 



The results for our model can be compared to the 'full diffusion' model, where 
nucleation and aggregation are due to the thermally activated diffusion process. 
Even if a detailed comparison is postponed to a future paper^^) which will extend 
our calculations and simulations to two dimensions and to a sequential diffusion 
model (see below), some comments are in order here. 

In Section [21 we already stressed an important difference concerning the small 
coverage Q behavior of island density, his ~ with an exponent x which is different 
for the two models (in d* = 2 the difference is even more relevant, because x = 2 
for our 'no diffusion' model, while x = 3 for the 'full diffusion' model). Another 
difference concerns the shape of the size distribution of islands. The origin of these 
differences should be traced back to Eq. (jl-lj) . which gives a rough evaluation of the 
probability that a third atom intervenes during the nucleation or aggregation process 
of a given atom. As a matter of fact, this criterion has two weak points, both related 
to the actual meaning of the diffusion length. 

First, in the early growing regime atoms may travel a distance larger than £d 
before being incorporated: can be correctly defined as a typical distance only in 
the regime where the island distance is approximately constant. Second, two atoms 
may stick together even if their 'initial' distance is larger then ij^, or, similarly, they 
may not stick even if their distance is smaller than £□: is an average quantity and 
taking it as a capture length kills diffusion noise. 

In order to evaluate the effect of the two features of our model, sequential- 
ity and absence of diffusion noise, we plan to study an intermediate model, which 
we can call 'sequential diffusion' model. This model is still sequential, but nucle- 
ation/aggregation does not occur deterministically. The capture length recipe is 
replaced by allowing a deposited atom to walk a fixed maximum number of ran- 
dom hops. Preliminary simulations^^^ in one dimension show that this model has 
statistical properties more similar to the 'full diffusion' model. 
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